The nucleoside analog 4'-fluorouridine suppresses the replication of multiple enteroviruses by targeting 3D polymerase.

Human enteroviruses are the major pathogens causing hand-foot-and-mouth disease in infants and young children throughout the world, and infection with enterovirus is also associated with severe complications, such as aseptic meningitis and myocarditis. However, there are no antiviral drugs available to treat enteroviruses infection at present. In this study, we found that 4'-fluorouridine (4'-FlU), a nucleoside analog with low cytotoxicity, exhibited broad-spectrum activity against infections of multiple enteroviruses with EC50 values at low micromolar levels, including coxsackievirus A10 (CV-A10), CV-A16, CV-A6, CV-A7, CV-B3, enterovirus A71 (EV-A71), EV-A89, EV-D68, and echovirus 6. With further investigation, the results indicated that 4'-FlU directly interacted with the RNA-dependent RNA polymerase of enterovirus, the 3D pol, and impaired the polymerase activity of 3D pol, hence inhibiting viral RNA synthesis and significantly suppressing viral replication. Our findings suggest that 4'-FlU could be promisingly developed as a broad-spectrum direct-acting antiviral agent for anti-enteroviruses therapy.

Predicting anti-SARS-CoV-2 activities of chemical compounds using machine learning models.

To accelerate the discovery of novel drug candidates for Coronavirus Disease 2019 (COVID-19) therapeutics, we reported a series of machine learning (ML)-based models to accurately predict the anti-SARS-CoV-2 activities of screening compounds. We explored 6 popular ML algorithms in combination with 15 molecular descriptors for molecular structures from 9 screening assays in the COVID-19 OpenData Portal hosted by NCATS. As a result, the models constructed by k-nearest neighbors (KNN) using the molecular descriptor GAFF+RDKit achieved the best overall performance with the highest average accuracy of 0.68 and relatively high average area under the receiver operating characteristic curve of 0.74, better than other ML algorithms. Meanwhile, The KNN model for all assays using GAFF+RDKit descriptor outperformed using other descriptors. The overall performance of our developed models was better than REDIAL-2020 (R). A web server (https://clickff.org/amberweb/covid-19-cp) was developed to enable users to predict anti-SARS-CoV-2 activities of arbitrary compounds using the COVID-19-CP (P) models. Besides the descriptor-based machine learning models, we also developed graph-based Attentive FP (A) models for the 9 assays. We found that the Attentive FP models achieved a comparable performance to that of COVID-19-CP and outperformed the REDIAL-2020 models. The consensus prediction utilizing both COVID-19-CP and Attentive FP can significantly boost the prediction accuracy as assessed by comparing its performance with other three individual models (R, P, A) utilizing the Wilcoxon signed-rank test, thus can ultimately improve the success rate of COVID-19 drug discovery.

In Silico Screening of Natural Flavonoids against 3-Chymotrypsin-like Protease of SARS-CoV-2 Using Machine Learning and Molecular Modeling.

The "Long-COVID syndrome" has posed significant challenges due to a lack of validated therapeutic options. We developed a novel multi-step virtual screening strategy to reliably identify inhibitors against 3-chymotrypsin-like protease of SARS-CoV-2 from abundant flavonoids, which represents a promising source of antiviral and immune-boosting nutrients. We identified 57 interacting residues as contributors to the protein-ligand binding pocket. Their energy interaction profiles constituted the input features for Machine Learning (ML) models. The consensus of 25 classifiers trained using various ML algorithms attained 93.9% accuracy and a 6.4% false-positive-rate. The consensus of 10 regression models for binding energy prediction also achieved a low root-mean-square error of 1.18 kcal/mol. We screened out 120 flavonoid hits first and retained 50 drug-like hits after predefined ADMET filtering to ensure bioavailability and safety profiles. Furthermore, molecular dynamics simulations prioritized nine bioactive flavonoids as promising anti-SARS-CoV-2 agents exhibiting both high structural stability (root-mean-square deviation < 5 Å for 218 ns) and low MM/PBSA binding free energy (<-6 kcal/mol). Among them, KB-2 (PubChem-CID, 14630497) and 9-O-Methylglyceofuran (PubChem-CID, 44257401) displayed excellent binding affinity and desirable pharmacokinetic capabilities. These compounds have great potential to serve as oral nutraceuticals with therapeutic and prophylactic properties as care strategies for patients with long-COVID syndrome.

ARPC5 acts as a potential prognostic biomarker that is associated with cell proliferation, migration and immune infiltrate in gliomas.

BACKGROUND: Gliomas are the most common malignant brain tumors, with powerful invasiveness and an undesirable prognosis. Actin related protein 2/3 complex subunit 5 (ARPC5) encodes a component of the Arp2/3 protein complex, which plays a significant role in regulating the actin cytoskeleton. However, the prognostic values and biological functions of ARPC5 in gliomas remain unclear. METHODS: Based on the TCGA, GEO, HPA, and UALCAN database, we determined the expression of ARPC5 in glioma. The results were verified by immunohistochemistry and Western blot analysis of glioma samples. Moreover, Kaplan-Meier curves, ROC curves, Cox regression analyses, and prognostic nomograms were used to observe the correlation between the ARPC5 expression and the prognosis of glioma patients. GO and KEGG enrichment analyses were conducted to identify immune-related pathways involved with the differential expression of ARPC5. Subsequently, the TCGA database was used to estimate the relationship between ARPC5 expression and immunity-related indexes, such as immune scores, infiltrating immune cells, and TMB. The TCIA database was used to assess the correlation between ARPC5 with immunotherapy. The association between ARPC5 and T cells marker CD3 was also evaluated through immunohistochemistry methods. The correlation between ARPC5 and T cell, as well as the prognosis of patients, was also evaluated using immunological methods. Moreover, the effect of ARPC5 on the biological characteristics of LN229 and U251 cells was determined by MTT, clone formation, and transwell migration assay. RESULTS: The high degree of ARPC5 was correlated with worse prognosis and unfavorable clinical characteristics of glioma patients. In the analysis of GO and KEGG, it is shown that ARPC5 was strongly correlated with multiple immune-related signaling pathways. The single-cell analysis revealed that ARPC5 expression was increased in astrocytes, monocytes and T cells. In addition, ARPC5 expression was strongly associated with immune scores, infiltrating immune cells, TMB, MSI, immune biomarkers, and immunotherapy. In experimental analysis, we found that ARPC5 was significantly overexpressed in gliomas and closely correlated with patient prognosis and CD3 expression. Functionally, the knockout of ARPC5 significantly reduced the proliferation and invasion of LN229 and U251 cells. CONCLUSIONS: Our study revealed that the high expression level of ARPC5 may serve as a promising prognostic biomarker and be associated with tumor immunity in glioma.

Discovery of Potent and Selective CB2 Agonists Utilizing a Function-Based Computational Screening Protocol.

Nowadays, the identification of agonists and antagonists represents a great challenge in computer-aided drug design. In this work, we developed a computational protocol enabling us to design/screen novel chemicals that are likely to serve as selective CB2 agonists. The principle of this protocol is that by calculating the ligand-residue interaction profile (LRIP) of a ligand binding to a specific target, the agonist-antagonist function of a compound is then able to be determined after statistical analysis and free energy calculations. This computational protocol was successfully applied in CB2 agonist development starting from a lead compound, and a success rate of 70% was achieved. The functions of the synthesized derivatives were determined by in vitro functional assays. Moreover, the identified potent CB2 agonists and antagonists strongly interact with the key residues identified using the already known potent CB2 agonists/antagonists. The analysis of the interaction profile of compound 6, a potent agonist, showed strong interactions with F2.61, I186, and F2.64, while compound 39, a potent antagonist, showed strong interactions with L17, W6.48, V6.51, and C7.42. Still, some residues including V3.32, T3.33, S7.39, F183, W5.43, and I3.29 are hotspots for both CB2 agonists and antagonists. More significantly, we identified three hotspot residues in the loop, including I186 for agonists, L17 for antagonists, and F183 for both. These hotspot residues are typically not considered in CB1/CB2 rational ligand design. In conclusion, LRIP is a useful concept in rationally designing a compound to possess a certain function.

Comparison and summary of in silico prediction tools for CYP450-mediated drug metabolism.

The cytochrome P450 (CYP450) enzyme system is responsible for the metabolism of more than two-thirds of xenobiotics. This review summarizes reports of a series of in silico tools for CYP450 enzyme-drug interaction predictions, including the prediction of sites of metabolism (SOM) of a drug and the identification of inhibitor/substrates for CYP subtypes. We also evaluated four prediction tools to identify CYP inhibitors utilizing 52 of the most frequently prescribed drugs. ADMET Predictor and CYPlebrity demonstrated the best performance. We hope that this review provides guidance for choosing appropriate enzyme prediction tools from a variety of in silico platforms to meet individual needs. Such predictions are useful for medicinal chemists to prioritize their designed compounds for further drug discovery.

Phosphorylation at Ser289 Enhances the Oligomerization of Tau Repeat R2.

In tauopathies such as Alzheimer's disease (AD), aberrant phosphorylation causes the dissociation of tau proteins from microtubules. The dissociated tau then aggregates into sequent forms from soluble oligomers to paired helical filaments and insoluble neurofibrillary tangles (NFTs). NFTs is a hallmark of AD, while oligomers are found to be the most toxic form of the tau aggregates. Therefore, understanding tau oligomerization with regard to abnormal phosphorylation is important for the therapeutic development of AD. In this study, we investigated the impact of phosphorylated Ser289, one of the 40 aberrant phosphorylation sites of full-length tau proteins, on monomeric and dimeric structures of tau repeat R2 peptides. We carried out intensive replica exchange molecular dynamics simulation with a total simulation time of up to 0.1 ms. Our result showed that the phosphorylation significantly affected the structures of both the monomer and the dimer. For the monomer, the phosphorylation enhanced ordered-disordered structural transition and intramolecular interaction, leading to more compactness of the phosphorylated R2 compared to the wild-type one. As to the dimer, the phosphorylation increased intermolecular interaction and ß-sheet formation, which can accelerate the oligomerization of R2 peptides. This result suggests that the phosphorylation at Ser289 is likely to promote tau aggregation. We also observed a phosphorylated Ser289-Na+-phosphorylated Ser289 bridge in the phosphorylated R2 dimer, suggesting an important role of cation ions in tau aggregation. Our findings suggest that phosphorylation at Ser289 should be taken into account in the inhibitor screening of tau oligomerization.

In silico binding affinity prediction for metabotropic glutamate receptors using both endpoint free energy methods and a machine learning-based scoring function.

Metabotropic glutamate receptors (mGluRs) play an important role in regulating glutamate signal pathways, which are involved in neuropathy and periphery homeostasis. mGluR4, which belongs to Group III mGluRs, is most widely distributed in the periphery among all the mGluRs. It has been proved that the regulation of this receptor is involved in diabetes, colorectal carcinoma and many other diseases. However, the application of structure-based drug design to identify small molecules to regulate the mGluR4 receptor is limited due to the absence of a resolved mGluR4 protein structure. In this work, we first built a homology model of mGluR4 based on a crystal structure of mGluR8, and then conducted hierarchical virtual screening (HVS) to identify possible active ligands for mGluR4. The HVS protocol consists of three hierarchical filters including Glide docking, molecular dynamic (MD) simulation and binding free energy calculation. We successfully prioritized active ligands of mGluR4 from a set of screening compounds using HVS. The predicted active ligands based on binding affinities can almost cover all the experiment-determined active ligands, with only one ligand missed. The correlation between the measured and predicted binding affinities is significantly improved for the MM-PB/GBSA-WSAS methods compared to the Glide docking method. More importantly, we have identified hotspots for ligand binding, and we found that SER157 and GLY158 tend to contribute to the selectivity of mGluR4 ligands, while ALA154 and ALA155 could account for the ligand selectivity to mGluR8. We also recognized other 5 key residues that are critical for ligand potency. The difference of the binding profiles between mGluR4 and mGluR8 can guide us to develop more potent and selective modulators. Moreover, we evaluated the performance of IPSF, a novel type of scoring function trained by a machine learning algorithm on residue-ligand interaction profiles, in guiding drug lead optimization. The cross-validation root-mean-square errors (RMSEs) are much smaller than those by the endpoint methods, and the correlation coefficients are comparable to the best endpoint methods for both mGluRs. Thus, machine learning-based IPSF can be applied to guide lead optimization, albeit the total number of actives/inactives are not big, a typical scenario in drug discovery projects.

Physiologically-Based Pharmacokinetics Modeling for Hydroxychloroquine as a Treatment for Malaria and Optimized Dosing Regimens for Different Populations.

Malaria is a severe parasite infectious disease with high fatality. As one of the approved treatments of this disease, hydroxychloroquine (HCQ) lacks clinical administration guidelines for patients with special health conditions and co-morbidities. This may result in improper dosing for different populations and lead them to suffer from severe side effects. One of the most important toxicities of HCQ overdose is cardiotoxicity. In this study, we built and validated a physiologically based pharmacokinetic modeling (PBPK) model for HCQ. With the full-PBPK model, we predicted the pharmacokinetic (PK) profile for malaria patients without other co-morbidities under the HCQ dosing regimen suggested by Food and Drug Administration (FDA) guidance. The PK profiles for different special populations were also predicted and compared to the normal population. Moreover, we proposed a series of adjusted dosing regimens for different populations with special health conditions and predicted the concentration-time (C-T) curve of the drug plasma concentration in these populations which include the pregnant population, elderly population, RA patients, and renal impairment populations. The recommended special population-dependent dosage regimens can maintain the similar drug levels observed in the virtual healthy population under the original dosing regimen provided by FDA. Last, we developed mathematic formulas for predicting dosage based on a patient's body measurements and two indexes of renal function (glomerular filtration rate and serum creatine level) for the pediatric and morbidly obese populations. Those formulas can facilitate personalized treatment of this disease. We hope to provide some advice to clinical practice when taking HCQ as a treatment for malaria patients with special health conditions or co-morbidities so that they will not suffer from severe side effects due to higher drug plasma concentration, especially cardiotoxicity.

In Silico Prediction of Pharmacokinetic Profile for Human Oral Drug Candidates Which Lack Clinical Pharmacokinetic Experiment Data.

BACKGROUNDS AND OBJECTIVES: In silico methods which can generate high-quality physiologically based pharmacokinetic (PBPK) models for arbitrary drug candidates are greatly needed to select developable drug candidates that escape drug attrition because of the poor pharmacokinetic profile. The purpose of this study is to develop a novel protocol to preliminarily predict the concentration profile of a target drug based on the PBPK model of a structurally similar template drug by combining two software platforms for PBPK modeling, the SimCYP simulator and ADMET Predictor. METHODS: The method was evaluated by utilizing 13 drug pairs from 18 drugs in the built-in database of the SimCYP software. All drug pairs have Tanimoto scores (TS) no less than 0.5. As each drug in a drug pair can serve as both target and template, 26 sets were studied in this work. Three versions (V1, V2 and V3) of models for the target drug were constructed by replacing the corresponding parameters of the template drug step by step with those predicted by ADMET Predictor for the target drug. V1 represents the replacement of molecular weight (MW), V2 includes the replacement of parameter MW, fraction unbound in plasma (fu), blood-to-plasma partition ratio (B/P), logarithm of the octanol-buffer partition coefficient (log Po:w) and acid dissociation constant (pKa). In V3, all above-mentioned parameters as well as human jejunum effective permeability (Peff), Vd and cytochrome P450 (CYP) metabolism parameters (Km, Vmax or CLint) are modified. Normalized root mean square error (NRMSE) was used for the evaluation of the model performance. RESULTS: We found that the performance of the three versions of the models depends on structural similarity of the drug pairs. For Group I drug pairs (TS ≤ 0.7), V2 and V3 performed better than V1 in terms of NRMSE; for Group II drug pairs (0.7 < TS ≤ 0.9), 8 out of 10 V3 models had NRMSE < 0.2, the cutoff we applied to judge whether the simulated concentration-time (C-T) curve was satisfactory or not. V3 outperformed the V1 and V2 versions. For the two drug pairs belonging to Group III (TS > 0.9), V2 outperformed V1 and V3, suggesting more unnecessary replacement can lower the performance of PBPK models. We also investigated how the prediction accuracy of ADMET Predictor as well as its collaboration with SimCYP influences the quality of PBPK models constructed using SimCYP. CONCLUSION: In conclusion, we generated practical guidance on applying two mainstream software packages, ADMET Predictor and SimCYP, to construct PBPK models for drugs or drug candidates that lack ADME parameters in model construction.

A multiple-step in silico screening protocol to identify allosteric inhibitors of Spike-hACE2 binding.

While the COVID-19 pandemic continues to worsen, effective medicines that target the life cycle of SARS-CoV-2 are still under development. As more highly infective and dangerous variants of the coronavirus emerge, the protective power of vaccines will decrease or vanish. Thus, the development of drugs, which are free of drug resistance is direly needed. The aim of this study is to identify allosteric binding modulators from a large compound library to inhibit the binding between the Spike protein of the SARS-CoV-2 virus and human angiotensin-converting enzyme 2 (hACE2). The binding of the Spike protein to hACE2 is the first step of the infection of host cells by the coronavirus. We first built a compound library containing 77 448 antiviral compounds. Molecular docking was then conducted to preliminarily screen compounds which can potently bind to the Spike protein at two allosteric binding sites. Next, molecular dynamics simulations were performed to accurately calculate the binding affinity between the spike protein and an identified compound from docking screening and to investigate whether the compound can interfere with the binding between the Spike protein and hACE2. We successfully identified two possible drug binding sites on the Spike protein and discovered a series of antiviral compounds which can weaken the interaction between the Spike protein and hACE2 receptor through conformational changes of the key Spike residues at the Spike-hACE2 binding interface induced by the binding of the ligand at the allosteric binding site. We also applied our screening protocol to another compound library which consists of 3407 compounds for which the inhibitory activities of Spike/hACE2 binding were measured. Encouragingly, in vitro data supports that the identified compounds can inhibit the Spike-ACE2 binding. Thus, we developed a promising computational protocol to discover allosteric inhibitors of the binding of the Spike protein of SARS-CoV-2 to the hACE2 receptor, and several promising allosteric modulators were discovered.

Aicardi-Goutières syndrome-associated mutation at ADAR1 gene locus activates innate immune response in mouse brain.

BACKGROUND: Aicardi-Goutières syndrome (AGS) is a severe infant or juvenile-onset autoimmune disease characterized by inflammatory encephalopathy with an elevated type 1 interferon-stimulated gene (ISG) expression signature in the brain. Mutations in seven different protein-coding genes, all linked to DNA/RNA metabolism or sensing, have been identified in AGS patients, but none of them has been demonstrated to activate the IFN pathway in the brain of an animal. The molecular mechanism of inflammatory encephalopathy in AGS has not been well defined. Adenosine Deaminase Acting on RNA 1 (ADAR1) is one of the AGS-associated genes. It carries out A-to-I RNA editing that converts adenosine to inosine at double-stranded RNA regions. Whether an AGS-associated mutation in ADAR1 activates the IFN pathway and causes autoimmune pathogenesis in the brain is yet to be determined. METHODS: Mutations in the ADAR1 gene found in AGS patients were introduced into the mouse genome via CRISPR/Cas9 technology. Molecular activities of the specific p.K999N mutation were investigated by measuring the RNA editing levels in brain mRNA substrates of ADAR1 through RNA sequencing analysis. IFN pathway activation in the brain was assessed by measuring ISG expression at the mRNA and protein level through real-time RT-PCR and Luminex assays, respectively. The locations in the brain and neural cell types that express ISGs were determined by RNA in situ hybridization (ISH). Potential AGS-related brain morphologic changes were assessed with immunohistological analysis. Von Kossa and Luxol Fast Blue staining was performed on brain tissue to assess calcification and myelin, respectively. RESULTS: Mice bearing the ADAR1 p.K999N were viable though smaller than wild type sibs. RNA sequencing analysis of neuron-specific RNA substrates revealed altered RNA editing activities of the mutant ADAR1 protein. Mutant mice exhibited dramatically elevated levels of multiple ISGs within the brain. RNA ISH of brain sections showed selective activation of ISG expression in neurons and microglia in a patchy pattern. ISG-15 mRNA was upregulated in ADAR1 mutant brain neurons whereas CXCL10 mRNA was elevated in adjacent astroglia. No calcification or gliosis was detected in the mutant brain. CONCLUSIONS: We demonstrated that an AGS-associated mutation in ADAR1, specifically the p.K999N mutation, activates the IFN pathway in the mouse brain. The ADAR1 p.K999N mutant mouse replicates aspects of the brain interferonopathy of AGS. Neurons and microglia express different ISGs. Basal ganglia calcification and leukodystrophy seen in AGS patients were not observed in K999N mutant mice, indicating that development of the full clinical phenotype may need an additional stimulus besides AGS mutations. This mutant mouse presents a robust tool for the investigation of AGS and neuroinflammatory diseases including the modeling of potential "second hits" that enable severe phenotypes of clinically variable diseases.

In silico binding profile characterization of SARS-CoV-2 spike protein and its mutants bound to human ACE2 receptor.

Severe acute respiratory syndrome coronavirus (SARS-CoV-2), a novel coronavirus, has brought an unprecedented pandemic to the world and affected over 64 million people. The virus infects human using its spike glycoprotein mediated by a crucial area, receptor-binding domain (RBD), to bind to the human ACE2 (hACE2) receptor. Mutations on RBD have been observed in different countries and classified into nine types: A435S, D364Y, G476S, N354D/D364Y, R408I, V341I, V367F, V483A and W436R. Employing molecular dynamics (MD) simulation, we investigated dynamics and structures of the complexes of the prototype and mutant types of SARS-CoV-2 spike RBDs and hACE2. We then probed binding free energies of the prototype and mutant types of RBD with hACE2 protein by using an end-point molecular mechanics Poisson Boltzmann surface area (MM-PBSA) method. According to the result of MM-PBSA binding free energy calculations, we found that V367F and N354D/D364Y mutant types showed enhanced binding affinities with hACE2 compared to the prototype. Our computational protocols were validated by the successful prediction of relative binding free energies between prototype and three mutants: N354D/D364Y, V367F and W436R. Thus, this study provides a reliable computational protocol to fast assess the existing and emerging RBD mutations. More importantly, the binding hotspots identified by using the molecular mechanics generalized Born surface area (MM-GBSA) free energy decomposition approach can guide the rational design of small molecule drugs or vaccines free of drug resistance, to interfere with or eradicate spike-hACE2 binding.

Drug-Drug Interaction Between Oxycodone and Diazepam by a Combined in Silico Pharmacokinetic and Pharmacodynamic Modeling Approach.

Opioids and benzodiazepines have complex drug-drug interactions (DDIs), which serve as an important source of adverse drug effects. In this work, we predicted the DDI between oxycodone (OXY) and diazepam (DZP) in the human body by applying in silico pharmacokinetic (PK) and pharmacodynamic (PD) modeling and simulation. First, we studied the PK interaction between OXY and DZP with a physiologically based pharmacokinetic (PBPK) model. Second, we applied molecular modeling techniques including molecular docking, molecular dynamics (MD) simulation, and the molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) free energy method to predict the PD-DDI between these two drugs. The PK interaction between OXY and DZP predicted by the PBPK model was not obvious. No significant interaction was observed between the two drugs at normal doses, though very high doses of DZP demonstrated a non-negligible inhibitory effect on OXY metabolism. On the contrary, the molecular modeling study shows that DZP has potential to compete with OXY at the same binding pocket of the active µ-opioid receptor (MOR) and κ-opioid receptor (KOR). MD simulation and MM-PBSA calculation results demonstrated that there is likely a synergetic effect between OXY and DZP binding to opioid receptors, as OXY is likely to target the active MOR while DZP selectively binds to the active KOR. Thus, pharmacokinetics contributes slightly to the DDI between OXY and DZP although an overdose of DZP has been brought to attention. Pharmacodynamics is likely to play a more important role than pharmacokinetics in revealing the mechanism of DDI between OXY and DZP.

Machine learning on ligand-residue interaction profiles to significantly improve binding affinity prediction.

Structure-based virtual screenings (SBVSs) play an important role in drug discovery projects. However, it is still a challenge to accurately predict the binding affinity of an arbitrary molecule binds to a drug target and prioritize top ligands from an SBVS. In this study, we developed a novel method, using ligand-residue interaction profiles (IPs) to construct machine learning (ML)-based prediction models, to significantly improve the screening performance in SBVSs. Such a kind of the prediction model is called an IP scoring function (IP-SF). We systematically investigated how to improve the performance of IP-SFs from many perspectives, including the sampling methods before interaction energy calculation and different ML algorithms. Using six drug targets with each having hundreds of known ligands, we conducted a critical evaluation on the developed IP-SFs. The IP-SFs employing a gradient boosting decision tree (GBDT) algorithm in conjunction with the MIN + GB simulation protocol achieved the best overall performance. Its scoring power, ranking power and screening power significantly outperformed the Glide SF. First, compared with Glide, the average values of mean absolute error and root mean square error of GBDT/MIN + GB decreased about 38 and 36%, respectively. Second, the mean values of squared correlation coefficient and predictive index increased about 225 and 73%, respectively. Third, more encouragingly, the average value of the areas under the curve of receiver operating characteristic for six targets by GBDT, 0.87, is significantly better than that by Glide, which is only 0.71. Thus, we expected IP-SFs to have broad and promising applications in SBVSs.

Incorporating structural similarity into a scoring function to enhance the prediction of binding affinities.

In this study, we developed a novel algorithm to improve the screening performance of an arbitrary docking scoring function by recalibrating the docking score of a query compound based on its structure similarity with a set of training compounds, while the extra computational cost is neglectable. Two popular docking methods, Glide and AutoDock Vina were adopted as the original scoring functions to be processed with our new algorithm and similar improvement performance was achieved. Predicted binding affinities were compared against experimental data from ChEMBL and DUD-E databases. 11 representative drug receptors from diverse drug target categories were applied to evaluate the hybrid scoring function. The effects of four different fingerprints (FP2, FP3, FP4, and MACCS) and the four different compound similarity effect (CSE) functions were explored. Encouragingly, the screening performance was significantly improved for all 11 drug targets especially when CSE = S4 (S is the Tanimoto structural similarity) and FP2 fingerprint were applied. The average predictive index (PI) values increased from 0.34 to 0.66 and 0.39 to 0.71 for the Glide and AutoDock vina scoring functions, respectively. To evaluate the performance of the calibration algorithm in drug lead identification, we also imposed an upper limit on the structural similarity to mimic the real scenario of screening diverse libraries for which query ligands are general-purpose screening compounds and they are not necessarily structurally similar to reference ligands. Encouragingly, we found our hybrid scoring function still outperformed the original docking scoring function. The hybrid scoring function was further evaluated using external datasets for two systems and we found the PI values increased from 0.24 to 0.46 and 0.14 to 0.42 for A2AR and CFX systems, respectively. In a conclusion, our calibration algorithm can significantly improve the virtual screening performance in both drug lead optimization and identification phases with neglectable computational cost.

How Well Does the Extended Linear Interaction Energy Method Perform in Accurate Binding Free Energy Calculations?

With continually increased computer power, molecular mechanics force field-based approaches, such as the endpoint methods of molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and molecular mechanics generalized Born surface area (MM-GBSA), have been routinely applied in both drug lead identification and optimization. However, the MM-PB/GBSA method is not as accurate as the pathway-based alchemical free energy methods, such as thermodynamic integration (TI) or free energy perturbation (FEP). Although the pathway-based methods are more rigorous in theory, they suffer from slow convergence and computational cost. Moreover, choosing adequate perturbation routes is also crucial for the pathway-based methods. Recently, we proposed a new method, coined extended linear interaction energy (ELIE) method, to overcome some disadvantages of the MM-PB/GBSA method to improve the accuracy of binding free energy calculation. In this work, we have systematically assessed this approach using in total 229 protein-ligand complexes for eight protein targets. Our results showed that ELIE performed much better than the molecular docking and MM-PBSA method in terms of root-mean-square error (RMSE), correlation coefficient (R), predictive index (PI), and Kendall's τ. The mean values of PI, R, and τ are 0.62, 0.58, and 0.44 for ELIE calculations. We also explored the impact of the length of simulation, ranging from 1 to 100 ns, on the performance of binding free energy calculation. In general, extending simulation length up to 25 ns could significantly improve the performance of ELIE, while longer molecular dynamics (MD) simulation does not always perform better than short MD simulation. Considering both the computational efficiency and achieved accuracy, ELIE is adequate in filling the gap between the efficient docking methods and computationally demanding alchemical free energy methods. Therefore, ELIE provides a practical solution for the routine ranking of compounds in lead optimization.

Entecavir as a P2X7R antagonist ameliorates platelet activation and thrombus formation.

Platelet activation is the primary cause of thrombosis. The P2X7 receptor (P2X7R) is a therapeutic target of thrombosis. However, it is still unknown whether P2X7R activation affects platelet thrombus. Our molecular docking results showed that entecavir as a P2X7R antagonist interacted perfectly with the human P2X7R (hP2X7R) in silico simulation studies. Furthermore, our experimental data revealed that entecavir could act as a P2X7R antagonist to exert cytoprotective effects against platelet activation via protecting mitochondrial function, improving lipid peroxidation and increasing antioxidant activity. Correlated with this, entecavir inhibited platelet aggregation, dense-granule secretion, P-selectin expression, integrin activation and Ca2+ increase. In experimental mouse model, entecavir could significantly inhibit arteriovenous thrombosis and prolong the bleeding time. Furthermore, we found that entecavir had no significant effect on prothrombin time (PT), activated partial thrombin time (APTT), thrombin time (TT), fibrinogen (FIB), mean platelet volume (MPV) and platelet counts (PLT). This study demonstrates that entecavir markedly prevents platelet activation and thrombosis through inhibiting P2X7R without affecting coagulation system. Therefore, entecavir may be a potential candidate for treating thrombosis disease.

Introducing Virtual Oligomerization Inhibition to Identify Potent Inhibitors of Aß Oligomerization.

Amyloid-ß (Aß) oligomers are known as the most toxic form of Aß peptides, and they are a major contributor to Alzheimer's disease. Therefore, developing antagonist screening methods for the formation of Aß oligomers is urgent and of great interest. In this study, we introduce virtual oligomerization inhibition (VOI), a novel virtual screening protocol that applies atomistic simulation to quantitatively investigate the ability of a ligand in interfering Aß oligomerization and the formation of Aß oligomers. Results from the VOI performance on six known inhibitors of Aß aggregation (brazilin, curcumin, EGCG, ELND005, resveratrol, and tacrine) are in excellent agreement with the results of expensive experiments. Moreover, VOI can reveal the mechanism and kinetics of the inhibition process at the atomistic level. VOI not only improves the efficiency of the antagonist screening for Aß oligomerization but also reduces the cost of performing the task. Attractively, the principle of VOI can also be applied to inhibitor screening for the aggregation of other amyloid proteins/peptides.

Prediction of the Binding Affinities and Selectivity for CB1 and CB2 Ligands Using Homology Modeling, Molecular Docking, Molecular Dynamics Simulations, and MM-PBSA Binding Free Energy Calculations.

Cannabinoids are a group of chemical compounds that have been used for thousands of years due to their psychoactive function and systemic physiological effects. There are at least two types of cannabinoid receptors, CB1 and CB2, which belong to the G protein-coupled receptor superfamily and can trigger different signaling pathways to exert their physiological functions. In this study, several representative agonists and antagonists of both CB1 and CB2 were systematically studied to predict their binding affinities and selectivity against both cannabinoid receptors using a set of hierarchical molecular modeling and simulation techniques, including homology modeling, molecular docking, molecular dynamics (MD) simulations and end point binding free energy calculations using the molecular mechanics/Poisson-Boltzmann surface area-WSAS (MM-PBSA-WSAS) method, and molecular mechanics/generalized Born surface area (MM-GBSA) free energy decomposition. Encouragingly, the calculated binding free energies correlated very well with the experimental values and the correlation coefficient square (R2), 0.60, was much higher than that of an efficient but less accurate docking scoring function (R2 = 0.37). The hotspot residues for CB1 and CB2 in both active and inactive conformations were identified via MM-GBSA free energy decomposition analysis. The comparisons of binding free energies, ligand-receptor interaction patterns, and hotspot residues among the four systems, namely, agonist-bound CB1, agonist-bound CB2, antagonist-bound CB1, and antagonist-bound CB2, enabled us to investigate and identify distinct binding features of these four systems, with which one can rationally design potent, selective, and function-specific modulators for the cannabinoid receptors.